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Abstract.  From  Hipparcos  proper  motions,  we  have  analyzed  the  lo¬ 
cal  kinematics  of  the  Milky  Way  from  the  young  0-B5  stars,  classical 
Cepheids,  and  the  late-type  K-M  giants.  Comparing  the  derived  results, 
the  kinematical  discrepancies  between  each  subset  of  stars  have  been  an¬ 
alyzed,  and  the  possible  reasons  for  the  kinematical  differences  have  been 
discussed. 


1.  Introduction 

In  the  early  part  of  this  century,  the  kinematics  of  stars  near  the  Sun  have 
been  known  to  provide  crucial  information  regarding  both  the  structure  and  the 
evolution  of  the  Galaxy.  Lindblad  (1927)  and  Oort  (1927)  developed  an  axisym- 
metric  model  of  the  rotation.  In  the  1950s,  the  stellar  kinematics  were  found 
to  vary  systematically  with  stellar  type.  Thereafter,  analyses  of  the  Galactic 
structure  and  kinematics  have  been  carried  out  via  various  observational  data. 
Considering  the  results  obtained  before,  IAU  has  recommended  a  set  of  Galactic 
constants  in  1985  (Kerr  &  Lynden-Bell,  1986).  Thanks  to  the  improvements 
of  the  observational  techniques,  many  important  results  for  Galactic  astronomy 
have  been  achieved  in  recent  years. 

The  present  work  will  concentrate  on  investigating  the  kinematical  features 
characterized  by  different  samples  and  discussing  the  possible  reasons  for  kine¬ 
matical  discrepancies  derived  from  the  different  stellar  samples. 


2.  Basic  Data  and  Analyses 

The  basic  materials  for  our  kinematical  analyses  of  the  Galaxy  include  astro- 
metric  parameters  of  stars  incorporated  with  astrophysical  measurements.  The 
Hipparcos  Catalogue  provides  observations  of  stars  with  precision  in  position, 
parallax,  and  proper  motion  of  1  mas  and  1  mas  yr-1. 

In  order  to  obtain  astrophysical  quantities  (MK  classification,  luminosity 
class,  spectroscopic  parallax,  etc.),  the  SKYMAP  catalog  containing  248,558 
stars  is  used.  Carrying  out  cross-identification  between  the  Hipparcos  and 
SKYMAP  catalogs,  excluding  all  binaries,  multiples,  and  suspected  non-single 
systems,  we  have  found  1,801  0-B5  stars  and  23,181  K-M  giants  common  to  both 


110 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

MAR  2000 

2.  REPORT  TYPE 

N/A 

3.  DATES  COVERED 

4.  TITLE  AND  SUBTITLE 

Stellar  Kinematics  from  Hipparcos  Proper  Motions 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Naval  Observatory  3450  Massachsuetts  Avenue,  N.W.  Washington, 
DC  20392-5420 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

Towards  models  and  constants  for  sub-microarcsecond  astrometry,  Proceedings  of  IAU  Colloquium  180 
held  at  the  U.S.  Naval  Observatory,  Washington,  DC,  USA,  27-30  March  2000 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

uu 

18.  NUMBER 
OF  PAGES 

5 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Stellar  Kinematics  from  Hipparcos 


111 


catalogs.  According  to  Westin’s  study,  the  Gould  belt  stars  are  embedded  in 
the  solar  vicinity  (Westin,  1985).  We  have  individually  removed  278  Gould  belt 
stars  from  our  0-B5  sample.  Hence,  the  sample  of  0-B5  stars  to  be  analyzed  is 
1,523  in  total. 

The  Hipparcos  Catalogue  provides  observations  of  223  classical  Cepheids. 
Based  on  the  spectroscopic  observations  (Metzger  et  al. ,  1998;  Pont  et  al. ,  1994), 
we  have  obtained  198  Hipparcos  Cepheids  with  radial  velocities.  In  the  present 
work,  we  adopted  the  PL  (period-luminosity)  relation  recently  proposed  by  Feast 
&  Catchpole  (1997)  to  calibrate  the  heliocentric  distances  of  Cepheids. 

The  first-order  Taylor  expansion  at  the  Sun  of  the  velocity  held  of  stars 
gives  a  three-dimensional  kinematical  model  (Clube,  1972).  The  equations  of 
condition  for  the  transverse  velocity  components  of  stars  contain  9  kinematical 
parameters  to  be  solved:  three  components  of  solar  motion  (uo,vo,wo),  three 
components  of  vorticity  (D32,  D21),  and  three  components  of  strain  velocity 

(H/ 2,  D is,  D32).  The  relevant  equations  were  given  by  Miyamoto  &  Soma  (1993). 

Combining  the  transverse  and  radial  velocities,  we  have  the  equations  of 
condition  with  12  parameters  to  be  determined,  including  three  additional  com¬ 
ponents  of  Galactic  expansion  and/or  contraction  (H//,  D22,  D33).  The  detailed 
description  is  given  by  Zhu  (2000).  Note  that  in  the  circular  stream  model  of 
Oort-Lindblad,  the  coefficients  Df2  and  D21  are  identical  to  the  familiar  Oort 
constants  A  and  B. 

In  a  galactocentric  cylindrical  coordinate  system  (R,9,z),  we  define  that  R 
is  measured  from  the  Galactic  center.  The  azimuthal  angle  9  is  reckoned  from 
the  a:  — axis  counterclockwise  around  the  axis  z.  Then,  we  have 


A  =  D+  = 
B  =  D21  = 
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where  Rq  (=8.5  kpc)  denotes  the  galactocentric  distance  of  the  Sun.  Vg  is  the 
rotational  velocity  of  the  local  standard  of  rest  (LSR)  about  the  Galactic  center. 

The  kinematical  parameters  will  be  solved  using  the  generalized  least-squares 
method.  Prior  to  applying  least-squares  to  the  equations  of  condition  to  derive 
the  parameters,  a  few  concerns  about  the  sampling  domain  of  stars  should  be 
raised. 

(1)  To  avoid  problems  at  the  solar  vicinity,  we  exclude  all  the  nearby  stars. 
For  the  0-B5  stars  and  Cepheids,  the  sampling  domain  is  confined  within 


|  z  |  <  0.35  kpc,  0.1  kpc  <  r  <  3.0  kpc.  (3) 

For  the  K-M  giants,  we  confine  ourselves  to  the  sample  within  the  domain 


z  |  <  0.5  kpc,  0.3  kpc  <  r  <  1.0  kpc.  (4) 


(2)  In  order  to  exclude  stars  with  extremely  erroneous  transverse  and  radial 
velocities  as  well  as  extremely  large  random  motions,  we  introduce  a  numerical 
filter  |  e  |.  Considering  both  the  measurement  errors  and  random  velocities,  we 
set  |  e  |=50  kpc  s-1  for  the  0-B5  stars  and  Cepheids,  and  |  e  |  =90  kpc  s-1  for 
the  K-M  giants. 
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3.  Solutions  of  Solar  Motion 

Solutions  of  solar  motion  for  each  type  of  stars  are  given  in  Table  1.  So  =  (uo,  Vo, 
wo)  is  the  total  velocity  of  the  Sun  relative  to  the  type  of  stars  with  the  apex 
towards  f@  and  b@. 


Table  1.  Solar  motion  derived  from  0-B5  stars,  classical  Cepheids, 
and  K-M  giants. 


Type  Data 

0-B5  PPMs 

Cep  PPMs 

Cep  PPMs+RVs 

gK-gM  PPMs 


uo  (km  s  1 ) 
11.59T0.49 
10.62±1.20 
9.79±0.86 
9.66±0.31 


v0 (km  s-1) 
13.39±0.48 
16.06il.14 

13.65i0.86 
21.45i0.32 


wo  (km  s  1 ) 
7.12i0.44 
8.60il.02 

8.40i0.86 
7.95i0.26 


So  (km  s 

19.09i0.48 

21.09il.17 

18.78i0.86 

24.87i0.30 


G(°)  bo(°) 

49.111.6  i21.9il.3 

56.0i3.5  i23.8i2.8 

54.4i2.9  i26.6i2.6 

65.810.7  il8.6i0.6 


The  solutions  of  solar  motion  derived  from  each  type  of  stars  are  apparently 
different.  The  radial  and  vertical  components  uo  and  wo  of  the  velocity  of  the 
Sun  do  not  vary  significantly  with  stellar  groups.  The  component  in  the  direction 
of  Galactic  rotation,  vo,  differs  from  type  to  type.  Possible  reasons  for  such  a 
large  variation  are  discussed  below. 

(1)  The  solar  motion  is  derived  from  individual  groups  of  stars.  It  is  well- 
known  that  the  solar  motion  with  respect  to  the  type  of  stars,  especially  the 
component  pointing  to  the  direction  of  Galactic  rotation,  is  strongly  correlated 
with  the  velocity  dispersion  (asymmetric  drift).  Analyzing  nearly  12,000  Hip- 
parcos  main-sequence  stars,  Dehnen  &  Binney  obtained  the  component  Vo  varied 
from  9  to  25  km  s-1  with  different  B-V  (1998).  The  youngest  disk  population  of 
stars,  the  0-B5  stars  and  Cepheids,  yields  much  lower  velocity  dispersion  than 
the  K-M  giants.  Thus,  the  solar  motion  Vo  derived  from  K-M  giants  should  be 
remarkably  larger  than  that  derived  from  0-B5  stars  and  classical  Cepheids. 

For  the  very  early  type  0-B5  stars  and  Cepheids,  the  variation  of  the  Vo 
component  cannot  be  explained  as  given  in  the  last  paragraph.  The  real  feature 
could  be  much  more  complicated. 

(2)  For  the  classical  Cepheids,  the  values  of  the  Vo  component,  derived  from 
Hipparcos  proper  motions,  from  radial  and  transverse  velocities,  are  significantly 
different.  The  causes  of  the  large  difference  in  solar  motion  from  Cepheids  should 
be  carefully  investigated  further.  One  of  the  possibilities  is  due  to  a  systematic 
shift  between  the  radial  velocity  system  and  the  proper  motion  system. 


4.  Galactic  Rotation  Solutions 

The  final  solutions  for  the  Galactic  kinematics  are  given  in  Table  2,  where  the 
unit  of  D  fj  is  km  s-1.  The  parameters  D^2  and  correspond  to  the  familiar 
Oort  constants  A  and  B  for  the  Oort-Lindblad  model,  given  in  eqs.  (1)  and  (2). 
Taking  Ro  =  8.5  kpc,  we  obtained  the  rotational  velocity  Vo,  listed  in  Table  2. 

The  IAU  has  recommended  a  set  of  Galactic  constants  in  1985,  in  which  Vo 
=  220  km  s-1  for  the  Galactic  rotation  and  Ro  =  8.5  kpc.  The  present  solution 
of  Galactic  rotation  are  remarkably  large,  compared  with  the  IAU  recommen¬ 
dation.  Analyzing  the  proper  motions  of  the  Southern  Proper  Motion  Program 
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Table  2.  Galactic  rotation  derived  from  0-B5  stars,  classical 
Cepheids,  and  K-M  giants.  The  unit  of  D±’s  js  km  s-1  kpc-1. 


Type  Data  D)"9  D91  D^,  D)"9  D99 

G-B5  PPiVIs  16.06il.13  -15.56±.83  1.89±.53  1.89±.53 

Cep  PPMs  16.08±1.14  -12.15i.86 

Cep  PPMsiRVs  15.37i0.85  -12.92i.85 

gK-gM  PPMs  15.51i0.93  -14.14i.75  -l.42i.90  -l.52i.64 


D+  Vo  (km  s  1 ) 

- 22 - 268.7ill.9 

240.0il2.1 
-2.60il.07  240.5il0.2 
252.1il0.1 


(SPM),  Mendez  et  al.  (1999)  derived  a  large  LSR  speed  of  270  km  s-1,  which 
is  in  excellent  agreement  with  ours  derived  from  the  0-B5  stars. 

The  difference  of  the  rotational  velocity  from  0-B5  stars  and  K-M  giants 
can  be  explained  with  the  velocity  dispersion  of  the  two  groups  of  stars.  The 
rotational  velocity  obtained  from  the  classical  Cepheids  is  obviously  slower  than 
that  of  K-M  giants  and  0-B5  stars.  One  of  the  main  causes  is  probably  due  to 
the  constituent  stars  of  moving  groups,  to  which  some  Cepheids  might  belong. 

A  meaningful  determination  of  parameters  D32  and  DJ2  (D32  =  DJ2  =  1.89 
±  0.53  km  s-1  kpc-1)  has  been  obtained  from  the  proper  motion  analysis  of  the 
0-B5  stars: 

{jw)r  =  ~2D 3"2  =  ~2D 3+2  =  _3'79  ±  1-05  km  S_1  kpc_1‘  (5) 


It  implies  a  clear  warping  motion  that  is  a  systematic  rotation  of  the  0-B5 
stars  about  the  axis  pointing  to  the  Galactic  center  in  the  sense  of  increasing 
the  inclination  of  the  H  I  warp,  which  we  have  discussed  in  our  previous  work 
(Miyamoto  &  Zhu,  1998).  Considering  the  solution  of  the  K-M  giants,  a  few 
points  will  be  added  here. 

The  K-M  giants  are  considered  to  be  well-relaxed  constituents  of  the  Galaxy, 
that  are  supposed  to  have  already  reached  a  steady  state  with  some  symmetry. 
The  present  solution  from  the  Hipparcos  K-M  giants  shows  that  the  K-M  giants 
seem  not  simply  to  exhibit  a  differential  rotation  along  the  Galactic  plane  as  one 
supposed,  because  of  a  meaningful  determination  of  the  parameters  D+g  and  D32. 
At  the  present  stage,  we  could  not  conclude  that  the  solution  of  the  parameters 
D?3  and  D32  reflects  a  new  class  of  Galactic  interior  motion  derived  from  the 
Hipparcos  proper  motions  of  the  K-M  giants.  If  we  suppose  the  determination 
of  Djg  and  DJ2  derived  from  the  K-M  giants  only  reflect  a  fictitious  rotation 
of  the  Hipparcos  proper  motion  system,  the  warping  motion  derived  from  the 
0-B5  stars  will  be  more  pronounced. 

The  solution  from  analyzing  the  proper  motions  and  radial  velocities  of  the 
classical  Cepheids  shows  that  Cepheids  seem  to  exhibit  a  contracting  motion 
in  the  solar  neighborhood  that  is  along  the  direction  pointing  to  the  Galactic 
rotation.  On  the  assumption  of  the  pure  axisymmetric  circular  rotation  of  the 
Galaxy  (V#  =  0),  we  have 


D 


+ 

22 


1  OVg 


Ro 


R  dd 


2.60  ±  1.07  km  s  1  kpc  1 . 


(6) 
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The  present  solution  seems  to  suggest  a  slight  contraction  of  the  Cepheids 
in  the  solar  neighborhood.  However,  the  real  feature  of  the  contraction  should 
be  carefully  investigated  further,  before  reaching  a  conclusive  explanation. 

One  of  the  possible  reasons  for  the  apparent  contraction  might  be  mem¬ 
bership  of  stars  in  moving  groups.  For  example  some  Cepheids  in  our  sample 
may  be  constituents  of  the  moving  groups.  Therefore,  the  kinematical  solution 
should  be  infected  by  such  constituent  stars  of  moving  groups,  and  might  lead 
to  a  systematic  deviation  in  the  determination  of  the  rotational  velocity  of  the 
Galaxy. 


5.  Conclusion 

The  results  of  the  solar  motion  relative  to  the  stellar  type  are  in  good  agreement 
with  those  given  by  Dehnen  &  Binney  for  the  main  sequence  stars.  We  deter¬ 
mined  very  stable  solutions  for  the  components  uo  and  wo-  The  component  Vo 
varies  in  a  wide  range  depending  on  type  of  the  stellar  age. 

The  Galactic  rotation  derived  from  the  0-B5  stars  (~  270  km  s-1)  is  re¬ 
markably  larger  than  the  IAU  recommended  value.  Analyzing  Hipparcos  proper 
motions  from  the  0-B5  stars,  we  have  identified  a  clear  stellar  warping  motion 
that  is  an  asystematic  rotation  dNz/dd/K  =  —  3.79  ±  1.05  km  s-1  kpc-1  of 
stars  around  the  axis  pointing  to  the  Galactic  center  in  the  sense  of  increas¬ 
ing  the  inclination  of  the  H  I  warp.  The  solution  from  the  classical  Cepheids 
seems  to  suggest  a  slight  contracting  motion  of  Cepheids  along  the  direction  of 
the  Galactic  rotation.  The  real  feature  of  the  contraction  should  be  carefully 
investigated  further. 
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